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ABSTRACT 

The infrared spectra of starburst galaxies are dominated by the low-excitation 
lines of [Nell] and [SIII], and the stellar populations deduced from these spectra 
appear to lack stars larger than about 35 M©. The only exceptions to this result 
until now were low metallicity dwarf galaxies. We report our analysis of the mid- 
infrared spectra obtained with IRS on Spitzer of the starburst galaxy NGC 1222 
(Mkn 603). NGC 1222 is a large spheroidal galaxy with a starburst nucleus that 
is a compact radio and infrared source, and its infrared emission is dominated by 
the [Nelll] line. This is the first starburst of solar or near-solar metallicity, known 
to us, which is dominated by the high-excitation lines and which is a likely host 
of high mass stars. We model the emission with several different assumptions as 
to the spatial distibution of the high- and low-excitation fines and find that the 
upper mass cutoff in this galaxy is 40 Mq-100 Mq. 

Subject headings: galaxies-starburst; galaxies: stellar content; galaxies: individual 
(NGC 1222, Mkn 603); infrared: galaxies 
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1. Introduction 

Star formation in starburst galaxies is often concentrated in compact, highly obscured 
regions. It is difficult to measure the stellar content of these sources. The stars cannot 
be studied directly, but only through the nebulae they excite. The proper use of nebular 
diagnostics has been controversial. In a stud y of all the galaxies with infrared spectra in the 



literature as of 2003, 



Rigbv fc Rieke 



(120041 ) concluded that the young stellar populations 
in galaxies of solar metallicity have relatively low T^s and that their IMFs, if Salpeterian, 
cannot extend to 100 Mq. They suggested that the upper mass cutoffs are < 40 Mq. Dwarf 
galaxies of much lower than solar metallicity may have higher mass stars; masses as high 
as 60 M© are needed to explain the spectra of the three such galaxies in their sample. This 
apparent lack of high mass stars in starbursts is a major puzzle of star formation studies. 
Is it real? Or is the apparent lack due to observational bias or problems in interpretation? 

We report here on the infrared spectra of the unusual SO galaxy NGC 1222 (Markarian 
603; D=34 Mpc (ho/71 Mpc)^^.) This galaxy is remarkable because it is only slightly, 
if at all, less metal-rich than our Galaxy, yet has the mid-infrared spectral signature of 
the most massive, hot stars, including a high ratio of IRAS 60-100/im flux. Our analysis 
of the galaxy is based on ground-based high-resolution maps of the radio and infrared 
continuum and near-infrared spectra, and on mid-infrared spectra obtained with the IRS 
Spectrometer on the Spitzer Space Telescope. We model the mid-infrared emission from 
this galaxy as arising in the superposition of a compact and bright nebula that contains 
the youngest stars and produces the high-excitation ionic lines, and a "cooler" spectrum 
region of more extended (^500 pc) emission. The line ratios are compared to the results 
of STARBURST99 stellar populations and MAPPINGS and CLOUDY photoionization 
models. The models are assumption-rich and the results are only limits. However, even 
with these caveats we can say that this galaxy has the highest upper mass cutoff yet known 
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for a galaxy that is not a metal-poor dwarf. 

In the next section we describe the observations and how we model the galaxy, in the 
Analysis section the STARBURST99 and photoionization grids, and in the Discussion, our 
conclusions. 



2. Observations 



2.1. Previous Work 



Petrosian fc Burenkov 



(119931 ) obtained optical spectra and images of the Markarian 
galaxy, NGC 1222. The galaxy has the optical appearance of a triple system. They found 
that the main source has a starburst nucleus of high surface brightness, superimposed on 
a fainter elliptical component of approximately 13 kpc by 9 kpc extent to the 25 mag 
arcsec"^ isophote. Within 2' (20 kpc) of the nucleus are two compact sources of high surface 
brightness, CI and C2, which they argue are dwarf galaxies interacting with the main galaxy. 
In spite of the proximity of the compact components, the overall isophotal morphology, 
even at large galactocentric radii, remains fairly symmetric. The optical spectra show that 



all three sources are photo-ionized, rather t 
the characteristic spectrum of a starburst. 



lan shock excited, and t hat th e main source has 



Petrosian fc Burenkov 



(I1993I ) found elemental 



abundances of O, Ne and N to be 12 + 



respectively. 



og[x/ H] 8.57 ± 0.09, 8.53 ± 0.026 and 7.35 ± 0.17 



Colziel. Doyon. fc Demerd (I2OOII ) obtained K-band spectra of NGC 1222 and 



report the Br 7 flux of the galaxy in 2 apertures: 23.5 ± 0.8 x 10~^^ergs~^ cm~^ in a large 
beam of 3"x9" and 16.9 ± 0.4 in a 3"x 3" beam, also consistent with a bright, compact 
starburst with an additional extended ionized component. Ha -|- [Nil] images of NGC 1222 



show that the ionized gas is 



( lArmus. Heckman. fc Miley 



brig 



itest at the locations of the compact continuum sources 



1990). 
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NGC 1222 is a strong infrared source and is in the Revised Bright Galaxy Sample 



( ISanders et al 



20031). Its IRAS fluxes are 0.59, 2.28, 13.06 and 15.41 Jy at 12, 25, 60 and 



lOOjUm respe c tively ; it has the high ratio of 60/100/im flux typical of active starbursts. 



Elfhag et al. 



(119961 ) searched for CO in NGC 1222 and did not detect any, so the molecular 



gas mass is n ot directly de t ected . However, dust emission in NGC 1222 has been detected 



at 850/i ni by 



VIA by 



Dunne et al 



Thomas et al. 



1^ usi 



20001 ) using SCUBA and it has been mapped in HI with the 



(120041 ) The logs of the atomic gas and dust masses are 9.38 and 



6.66; NGC 1222 is at the low end of the sample for both quantities. 



2.2. Radio Continuum Maps 



The spatial distribution of the starburst is a key element of the interpretation of 
the mid-infrared spectra, and NRAO VLAl^ maps allow us to map the structure of the 
starburst. We show radio maps at 6 cm (C band) and 1.3 cm (K band) obtained from the 
VLA archives in Figures la and lb. The 1.3 cm image is shown superimposed on the 6 cm 
in Figure Ic. Details of the observations are in Table 1. The beam sizes of the observations 
differ by about a factor of two; for the overlay in Figure Ic, the maps have been convolved 
to the same 1.75" x 1.2" beam. The VLA maps are sensitive to structures less than ~15" in 
size. 



In the 6 cm map, there is a secondary source about 10" east of the main source, and 
there is a 3a feature in K-band at th e same location. Comparing the radio data to the 



only image of 



Petrosian fc Burenkov 



(119931 ) ■ we associate the second radio source with 



their CI object, which they identify as a separate interacting system. We see no radio 



^The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation, operated under cooperative agreement by Associated Universities, Inc. 
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counterpart to their C2 to the southwest of th e main source, a . 



[Nil] source (lArmus. Heckman. fc Miley 



199C 



Condon et al 



thoug h it is a bright Ha + 



fll990h find a size of 11" x8 



" at 20 cm and gives a total 20 cm fiux of 55 mJy in 18" and 39.7 in 5". The CI source 
m ay appear as an extens i on in their map. The position of our radio source agrees with that 



of 



Sramek fc Weedman 



( 1l986l ) 20 cm map and is about 10" off the quoted optical source 



position, but that position is quoted with large (±5") uncertainties. 

The radio emission at 1.3 cm comes almost entirely from a bright source about 4" (620 
pc) in diameter; the 1.3 cm emission has an unusual shape which could be an incomplete 
ring, or two bright peaks not completely separated. The 6 cm emission has the same 
bright source with structure consistent with the 1.3 cm shape; the beam size at 6 cm is 
almost twice that at 1.3 cm so the same detail cannot be seen. The 6 cm map also shows 
extended emission, mostly in the east-west direction. The radio maps are consistent with 
the 11.7/im image in the Keck/LWS images of Turner, Gorjian & Beck 2007, which also 
shows a bright, concentrated source extended EW. Neither CI nor C2 fall in the field of the 
LWS observations. 

Comparison of the VLA maps at 6 and 1.3 cm, and with single dish fiuxes, reveals 
much about the spatial distribution of the radio emission and the starburst. The galaxy 
looks very different in the two images, so there are clearly variations in spectral index 
and sources of the radio emission. The expected sources of radio emission in starburst 
galaxies are nonthermal synchrotron emission, with spectral index, a ~ —0.7 to —0.8, and 
thermal free-free emission from HII regions, with a ~ —0.1 [S oc z/"). The 1.3 cm source 
has flux, found using the AIPS program TVSTAT, of S\,j,cm = 14 ± 2 mJy in 28 square 
arcseconds. The central source at 6 cm, excluding the extended emission components, has 
'S'e cm = 14 ± 4 mJy in the same region. The central source is thus a bright 1.3 cm emitter, 
roughly as bright at 1.3 cm as at 6 cm, and the radio emission in the central region is 
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mostly if not entirely thermal. Single dish fluxes measured by iMarx et al. I (119941 ) with the 
Effelsberg 100 m telescope are 30.7 ±3 mJy at 6 cm and 14.8 mJy ±3 at 2.8 cm; the 100 m 
beamsizes are 145" at 6 cm and 69" at 2.8 cm. Comparison of the total mapped flux in 
our 6 cm VLA map to the Effelberg fluxes indicates that there is a significant extended 
component to the 6 cm emission that is not detected in our VLA maps. This extended 
component appears to be dominated by nonthermal synchrotron emission, based on the 
single dish spectral index of -0.9. By contrast, the 1.3 cm VLA flux is essentially equal to 
the 2.8 cm Effelsberg flux. Since optically thin thermal free-free emission from Hu regions 
has a nearly flat spectrum, with a z/""'^ spectrum, then the VLA map must detect nearly 
all of the 1.3 cm emission in this galaxy, and the source must be largely free- free emission. 
Although a strong thermal radio source, this central source is not t he brightest source in 



the image of Ha + [Nil] image of lArmus. Heckman. fc Mileyl (jl990l ). probably because of 
extinction (§3.2). 

Thus the high frequency radio fluxes indicate that the 1.3 cm VLA map of NGC 1222 
is mostly thermal free-free emission from HII regions, and that this emission is confined 
to the source shown in the map of Figure lb. In the VLA maps, the 1.3 cm emission in 
NGC 1222 appears even stronger relative to the 6 cm emission than the slight negative 
spectral decrement of thermal emission would require. This source may be one of those, 
common in extreme starbursts, where the radio spectrum rises inste ad of falling at high 
frequencies because lower frequencies have significant optical depth (ITurner. Ho. fc Beck 



1998 



Beck. Turner, fc Kovo 



2000 : 



Johnson et al. 



20011). 



The inner region of NGC 1222 contains a large population of young ionizing stars. If 
we assume the 1.3 cm emission in the VLA maps is entirely thermal free-free emission, 
the observed flux implies Nljc of 1.45 x 10^^ s^^ for the central 3". This is the ionization 
equivalent of 1.45 x 10^ standard 07 stars. 
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2.3. Spitzer IRS Spectra 

NGC 1222 was observed with the IRS on Spitzer in staring mode as part of the IRS 
Standard Spectra program, Program 14, for which J. Houck was Principal Investigator . 
The middle-infrared spectral diagnostics we need for this study fall in the bandpass of the 
Short-Hi module, which covers the wavelength range 9.9-19.6/im with spectral resolution of 
~ 600. The NGC 1222 data were obtained in two positions separated by about 3" roughly 
parallel to the slit. Their full Spitzer identifiers are in the Table; here they will be called 
Spectrum 1 and Spectrum 2. Spectrum 1 is shown in Figure 2. The line fluxes obtained for 
the diagnostic lines, for each spectrum, are shown in Table 2. We extracted the spectra from 
the post-BCD data provided by the Spitzer pipeline. IRS data requires different treatment 
for extended and point sources. Since NGC 1222 is n either a poi nt source nor a flat and 



20061 ). who use the extended 



smooth extended source, we follow the SINGS group (iSeth 
source option in the SPICE package to obtain line over continuum fluxes in clumpy galaxies. 
They have shown that where the flux calibration can be constrained by other data, the 
extended source calibration is usually accurate to within 20%. The uncertainty in the fluxes 
and the line ratios is dominated by the calibration difficulties: since the correction for 
spatial structure depends strongly on the wavelength, both the relative and the absolute 
calibration are affected. We assume 20% as the working uncertainty for all the lines. 

The nomimal position of Spectrum 2 is 03^ 08™56.7*, — 2°57'18" and of Spectrum 1 
03^ Qgm 55.9^^ —2° 57' 20"; the former agrees better with the radio position. However, 
the beam size in the IRS is large, on these scales, and also depends on wavelength: the 
beam is around 3" at 10/im and 6" at 20/im, so the positions actually overlap for most 
of the wavelength range. Spectrum 2 has somewhat stronger lines than Spectrum 1 but 
the difference is less than the calibration uncertainty in all cases, and the line ratios on 
which the analysis in the next section depends differ by less than 10% between the two 
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observations. We have therefore taken the average of the hne ratios for each observation 
and work with that number for simphcity. 



3. Analysis 



3.1. Spatial Stucture 

The ground-based radio, K-band and 11.7 fim continuum data give a consistent picture 
of NGC 1222: the young stars are concentrated into a central source no larger than 3" 



diameter, ~ 470pc, which produces most of the in 
Th is is the source called t h e star burst nucleus by 



rared and radio emission in the galaxy. 



by 



Petrosian fc Burenkov 



Colziel. Doyon. fc Demerd (120011 ) and 



(119931 ). We do not know if it is the kinematic nucleus of the 
galaxy, as there is little dynamic information, but it is clearly the dominant region of star 
formation and the actual starburst. (Note that the Spitzer Short-Hi observations do not 
cover the optical and 6 cm source CI, and we cannot judge whether it is also an RISN or 
some other kind of emission region.) There are more young stars in a less dense region 2-3 
times as large as the compact source. The phenomenon of a very concentrated group of 
young stars with strong radio and dominant infrared emission is a common one in starburst 
galaxies: w'e haye called such sources in other galaxies RISN (Radio-Infrared Super-Nebulae 



Beck et al. 



20021 ). 



How are we to analyse the infrared spectra in light of the observed spatial structure of 
the galaxy? The suite of infrared diagnostic lines can be for convenience thought of as the 
high-excitation lines, [Nelll] and [SIV], and the low-excitation lines [Nell] and [SHI]. (The 
[Aril] and [Arlll] fine-structure lines at 6.99 and 8.99 /im are also useful diagnostics, but 
do not fall in the wavelength region measured by the Spitzer high resolution modules, and 
cannot be seen clearly in the low- resolution results). The homo- and hetero-nuclear line 
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ratios give a measure of the relative strength in the total radiation field of the hard photons 
that create the ions of A^e++ and 5*+^ and the softer photons responsible for the other ions. 
From this measure of the ionization it is possible to work backwards to stellar populations 



that could produce that ionizing field. This is an active area of research and t 



nere are 



several programs and approaches, comprehensively reviewed in 



Rigby fc Rieke 



(120041 ) 



For the line ratios to be meaningful they must compare lines from the same source. 
The galaxy appears in both the radio and optical maps to have a compact main source 
embedded in more diffuse emission. The IRS slit is large enough to encompass most of the 
total emission, compact and diffuse together (Figure 1). If the high- and low-excitation 
lines are distributed differently between the compact and diffuse components, the global 
line ratios will be subject to misinterpretation. This phenomenon has been observed in 
other galaxies, particu larly those that are forming super star clusters, such as NGC 5253 



( ICrowther et al 



19991 ). It is also what we should expect, from the nature of the IMF. 



In any IMF the most massive stars will be only a small fraction of the total by number. 
To have enough very massive stars that they significantly affect the spectrum, the total 
number of stars must be very large. So the most massive stars will be found in NGC 1222, 
if at all, in the compact source, which as shown above must contain tens of thousands of O 
stars. A star must be larger than 40-45 Mq to excite as much [Nelll] as [Nell] emission. 
[Nell], in contrast, is excited by the more common late-0 and early B stars: even a BO star 
of only 20Mq can produce strong [Nell] emission. These lower mass stars are a much larger 
fraction of the total by number and will be common even in the sparser stellar population 
of the diffuse source. 

So we argue that almost all the [Nelll] should be assigned to the compact 3" source 
visible in the VLA 1.3 cm image. But what about the [Nell] and the sulfur lines? What 
fraction of the line strengths should be attributed to the compact source and what 
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to the extended? We note that IColzieL Doyon. &: Demera (l200ll ) found an excess of 
6.6 X lO^^^ergs"^ cm"^ of Br 7 emission in their large beam over that in their small. Their 
small beam being comparable in size to the compact source, and their larger beam to the 
entire Spitzer slit, we take this extended component of emission and estimate the [Nell] 
flux that should be produced by^ t hat (j uantity of ionized gas. We use the neon abundance 
found by jpetrosian fc Burenkov! (ll993i ) of 3.34 x 10~^ relative to H and assume that all 
the neon is singly ionized; the latter is the default case in normal H ll regions. With these 
assumptions the extended component is expected to produce a [Nell] line of 1.8 Jy strength. 
Which equals, within the calculational uncertainties, the entire [Nell] line flux seen by 
Spitzer! This argues that the low-excitation line of [Nell] is almost entirely produced in the 
extended emission and that the contribution of the compact source is negligible. 



So NGC 1222 appears to resemble NGC 5253 (ICrowther et al.lll999f ). where all the 
[Nelll] flux seen in the much larger beam of ISO proved to be emitted by a single compact 
source, and all the [Nell] from an extended emission region. How will the other lines track? 
We don't know enough of the structure of the ionization field to treat this fully, but can 
get a crude and approximate answer from the ionization potentials of the sulfur ions. S^~^ 
can only exist in the A^e^ zone, but S^^ can coexist with A^e"*" as well as with Ne'^^. So 
the extended diffuse emission could produce some part of the [SIV] as well as essentially all 
of the [SHI]. In NGC 5253, this was not the case, and all the [SIV] was produced by the 
compact source, but in II Zw 40 (another lo w-metallicity dwarf) as much as 3 0% of the SIV 



may come from outside the compact source (IMartin-Hernandez et al. 



20061). 



Since we do not know the spatial distribution of the mid-infrared lines, we therefore 
work with the line ratios we would derive for the following three cases. These cases should 
cover all possibilities for the distribution of the mid-infrared lines in NGC 1222. 



Model 1: no spatial separation; the lines are co-extensive in origin. We know from the 
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Br'y results and all the arguments above that this model is unrealistic. Nevertheless, we 
present it as an extreme lower limit to the stellar temperatures and masses. 

Model 2: the most stratified model; the high and low-excitation lines are essentially 
disjoint. In this case NGC 1222 would resemble NGC 5253, where all the low-excitation 
lines were from the diffuse component. For the calculations we assume that 90% of the 
[Nellll] and [SIV] are from the compact source and 90% of the [Nell] and SIII] from the 
diffuse. We chose 90% because the remainder it leaves is consistent with the weakest Spitzer 
lines and the calculated [Nell] flux. This will give the hottest and most massive stars. It is 
also, we think, the most realistic model. 



Model 3: an intermediate c ase, found by applying t 



le solar neon abundance rather 



than the higher value found by 



Petrosian &: Burenkov 



(I1993I ) to the gas producing the 



extended Brj emission. This gives 0.96 Jy of [Nell] in the extended emission, for a model 
in which 60% of the low excitation lines [SIII] and [Nell] are in the extended source and 
40% in the compact. In this model as in model 2, 90% of the high-excitation lines are from 
the compact source; the distributions of [Nell] and [SIII] are changed. 

These physical models are simple but they are consistent with everything we know 
about the source and about similar star formation regions and they will show how the 
stellar types deduced depend on the spatial distribution. Line ratios for the models are 
given in Tables 3 and 4. 



3.2. Extinction Effects 



In the above discussion we used 
extinction. Based on optical spectra. 



;he observed line strengths w ithout correction for 



Petrosian &: Burenkov 



I993I ) found the maximum 



reddening to be 0.4 mag, for a visual extinction of 1.2 mag. The i?r7/radio ratio tells 
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a different story; the Br^y ffuxes are considerably lower than the value predicted from 
the radio continuum flux, which for Tg = 7500 K and no extinction at Br'-f would give 
1.6 X 10^^'^erg s~^ cm~^ in 3"x3". The discrepancy with the observed gives 2.4 mag of 
extinction at B7 and 19 magnitudes Ay. The different extinctions derived from different 
markers are common in star formation regions, where the optical emission region and the 
infrared sources can be physically distinct, and where much of the extinction is produced 
within the HII regions themselves. 

An obscuration of 19 magnitudes is typical of very dense compact Hjj regions. If we 

1 have to 



Giveon. Morrisset. fc Sternberg 



(12002) and 



accept that value for the region producing the infrared lines in NGC 1222 we w i 
consider the effect on our line ratios. We follow 
take Asiv = 0.7Ak, A^eii = 0.28^4^, AncIii = 0.2274^, and Asm = 0.32^4^. The extinctions 
at the [Nell], [Nelll], and [SIII] lines are so similar that even with this high figure for the 
total extinction the line ratios will change by only about 20%. The uncertainty introduced 
by the extinction is thus less that that due to the spatial structure, the metal content, and 
the instrument calibration. The problem is more severe for [SIV], where the line fluxes 
would be increased by a factor of 4.8. Correcting the [SIV]/[SIII] and [SIV]/[NeII] ratios for 
this extinction would increase them by factors of 2.3 and 2.5, respectively. This extinction 
factor would move the line ratios calculated for model 1 close to those of model 3, and those 
of model 3 close to model 2. 

The distribution of obscuring matter in NGC 1222 is probably not uniform, just as the 
emission is not uniform. The extinction may be higher towards the compact source and 
lower in the extended emission region. Until the galaxy is mapped with high resolution in 
the infrared, the above results are only approximate. 
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3.3. Photoionization Models 



Deducing the stellar population from the observed infrared lines is a complicated 
process. The line ratios are compared to those predicted by models, and the inputs into 
the models, which include the time dependance of the star formation process (i.e., is it an 
instantaneous burst or continuous), the mass function and mass limits of the burst, the 
stellar atmospheres, and the metal content of the stars and of the nebula. Much work has 
gone into the establishment and testing of the photoionization models and improvements 
continue, especially on the stellar atmospheres, the inclusion of stellar winds, and the 
inclusion of the Wolf-Rayet stage. The models are constantly changing and improving. 



We use the STARBURST99 software and data package (ILeitherer et al 



19991 ) to 



model the ionizing spectra produced by different candidate stellar populations, and the 



MAPPINGS code flKewlev et al. 



results general 



by CLOUDY flFerland 



y in good agre ement ( 



20071) to calculate the out put spectrum. This method gives 



Rigbv fc Rieke 



20041 ) with output spectra calculated 



20061 ) from ionizing spectra generated by STARBURST99. The 



models are most limited, or least realistic, in their geometrical simplicity — real star 
formation regions are not likely to be spherical, symmetric or uniformly filled. In NGC 
1222, as in almost all galaxies, we know that there is complex structure but it can't be 
observed directly. So we use the default geometrical assumptions in the photoionization 
code of sphericity, uniform filling, and isobaric density structure with P/k = 10^ and mean 
temperature 10^. The ionization parameter is set by the STARBURST99 output luminosity 
and the density. The models compare bursts of 10^ Mq with a Kroupa IMF with power-law 
exponents 1.3 for 0.1 to 0.5 M© and 2.3 for higher masses. We assume Pauldrach/Hillier 
stellar atmospheres, Geneva evolution tracks with high mass-loss rates, and dust of standard 
depletions. We vary the parameters 1) upper mass limit of the IMF, which may be refered 
to as the mass and 2) the age. We also ran models for different values of metal content. 
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The infrared diagnostic lines are very sensitive to the metal content of the source. 
Starbursts with lower metal content will have hotter nebulae than others of the same mass 
parameters and ages but with more metals. This is partly because stars with fewer metals, 
and thus less line-blanketing, will have hotter spectra than metal rich stars of the same 
mass, and partly because the ionized nebula cool by emitting metal lines. Until now, the 
starbursts in which [Nelll] is clearly stronger than [Nell] have been low metal-content blue 
dwarf galaxies: NGC 5253, II Zw 40, and NGC 55. The m etal content of NGC 1222, in 



contrast, is close to Galactic: 



Petrosian &: Burenkov 



(I1993I ) describe it as "only slightly, if 



at all deficient", with [0/H] 70% of solar, [Ne/H] almost twice solar, and only [N/H], at 
35% of solar, really deficient. Other elements were not measured. STARBURST99 offers 
stellar atmospheres in a range of metallicities, but the only realistic options for this galaxy 
were 0.4 solar, solar or twice solar. We used solar value as it is closest to the reported. The 
metal content of the nebula can be adjusted in the codes and we used gas with the stated 
abundances of O, Ne and N and the solar abundance of all others. 

A weakness of the models is that the stellar atmosphere abundance cannot be 
fine-tuned, and if the stars in NGC 1222 have metal content significantly lower than solar, 
they will produce more high-excitation gas than predicted and the stellar types derived 
from the models may be misleadingly high. We tested the possible impact of this by 
running several models with stellar atmospheres of 0.4 solar metallicity. If the nebular gas 
abundance patterns were the same, the two stellar atmosphere models produced line ratios 
within 15% of each other. So unless the stellar metal abundances are radically different 
from the estimated, this will not strongly affect the results. 
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4. The Most Massive Stairs in the Stairburst 

4.1. The Compact Source 

In Figure 3 we show the inferred hne ratios we would observe for the three models 
of spatial distribution, overplotted on simulation results. For Model 1, which treats the 
whole galaxy as one uniform source of all mid-infrared lines, and will give the coolest, least 
massive stars, the [NeIII]/[NeII] ratio can be fit by age less than 2.2 Myr and an upper 
mass limit of 35 M©, or alternatively, age 2.2-3.2 Myr and any mass above 40 M©. In 
Model 2, which we favor because of the radio morphology, the extreme concentration model 
in which nearly all the low excitation emission comes from extended (>500 pc) gas, only 
an upper hmit of > 50 — 70 Mq stars and age less than 2 Myr can match the results. For 
Model 3, the intermediate case in which 40% of the low excitation [SIII] and [Nell] flux is 
emitted by the nuclear region, the neon line ratio gives an upper mass limit around 40Mq 
for age less than 2 Myr or up to 100 Mq if in the 2-3 Myr range where the ionization is 
coolest. 

[NeIII]/[NeII] has two clear advantages over the other hne ratios we can form with the 
Spitzer data: it is homonuclear, and the two lines have very similar, very low, extinction 
(the line ratio does not have to be corrected for extinction at the current level of analysis). 
The [SIV] line is very close to the silicate absorption feature and has the highest extinction 
of the useful infrared lines, and the sulfur line ratio, since it covers the greatest wavelength 
range, is also the most susceptible to the wavelength dependance of the extended source 
calibration. With these warnings, the [SIV] /[SIII] line ratio without extinction correction 
is best fit by, in Model 1: masses less than 30 Mq and ages less than 3 Myr; Model 2: 
masses greater than 50 Mq and ages less than 2.5 Myr, and Model 3: masses of 30 Mq and 
ages less than 3 Myr. Extinction corrections move these results towards higher masses; the 
corrections in section 3.2 will give mass limits of about 30, 40, and >70 Mq for the three 
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spatial models. 

The [SIV]/[NeII] line ratio can be affected both by the significant extinction and by 
abundance. The simulations show that model 1 can be fit with 30-40 M© stars less than 
2 Myr old, model 2 by >50 M©, and model 3 by a mass >40 M©. The extinction correction 
results moves the best fit to 40 and 70 Mq for models 1 and 3 and more than 100 Mq for 
model 2. 

Note that in many cases the line ratios can agree with either a very early age, or a time 
later than about 3.5 Myr when stars have entered the Wolf-Rayet phase, which produces a 
hard ionization field. We did not consider the second stage because the Wolf-Rayet feature 
has not been seen in NGC 1222. But the high extinction and the paucity of deep optical 
spectra in this galaxy mean that Wolf-Rayet stars have not been strictly ruled out, and the 
possibility should be kept in mind. 

4.2. The Diffuse Emission 

The [NeIII]/[NeII] ratio in the diffuse region, although based on assumptions, is 
plausible because the [Nelll] is almost certain to be strongly concentrated and because 
the [Nell] fiux is derived in a fairly straightforward way. The [SIII] cannot be so easily 
calculated because the abundance of sulfur is not known. With these considerations, the 
line ratios agree with a relatively low upper mass limit of 30 Mq, and age less than 4 Myr. 
In this regime. Models 2 and 3 fit the same photoionization grids equally well. 
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5. Discussion 

Deducing the stellar population of a star cluster from the nebular line it excites is 
never straightforward; we must find the best values of mass and age with only three line 
ratios observed. In NGC 1222 it is even more complex, because of the large slit on Spitzer 
which combined flux from the compact source and from the extended emission, so that we 
must dis-entangle the spatial distribution from the observed line ratios. We have presented 
many models and must now consider which are the most physical and probable. 

That the compact source in NGC 1222 must contain most of the [Nelll] and [SIV], 
and that the extended Br 7 emission must be associated with some [Nell] and [SIII], is 
certain. Reality is thus somewhere between our Models 2 and 3 (the most concentrated and 
intermediate cases). Among the line ratios, we give the highest weight to [NeIII]/[NeII] as 
it is homonuclear and not strongly affected by extinction. (The [SIV] /[SIII] ratio compares 
two lines whose critical densities are significantly different, 3.7 x 10^ for [SIV] and 9 x 10^ 
for [SIII], both within the range of densities expected for compact star formation regions, 
in addition to the considerable extinction correction.) So the best-fitting models are those 
with upper mass limit between 40 Mq and 100 Mq and age less than 2 Myr. These results 
cannot be better refined until high resolution observations of the infrared lines, which would 
directly measure the contributions of the compact and diffuse component, can be carried 
out. 



5.1. Comparision to Other Simulation Results 



We compare our model results to the CLOUDY simulations of iRigby fc Rieke 



(120041 ). 



They used the same STARBURST99 input spectra but CLOUDY instead of the Mappings 
photoionization code; in Figure 11 of their paper it may be seen that for high masses and 
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solar metallicity they agree very well. We show in Figure 4 the line ratios as a function 
of time for CLOUDY simulations and solar metallicity, taken from that paper. The 
[NeIII]/[NeII] line ratio of NGC 1222 is so high that even the unrealistic, least concentrated 
Model 1 would require an upper mass limit of ~ 50 Mq, and for the other models even 
100 M0 would not provide enough ionization. The [SIV]/[SIII] and [SIV]/[NeII] ratios 
behave similarly. 



Rigby fc Rieke 



he MAPPINGS models of this paper and the CLOUDY models of 
(I2OO4J ) used different parameters and structures for the ionized regions, but it is apparent 
that the nebular metallicity is the driving factor that makes the results so different. If the 
metallicity of NGC 1222 were 0.2 solar, it would be like NGC 5253, which has similar line 
ratios, and which (Figure 4 of Rigby & Rieke) can be fit with Mupper > 4OM0 . In a solar 
metallicity system, the line ratios of NGC 1222 can only be fit by truly extreme stellar 
populations. In the moderately sub-solar metal models which are most appropriate for the 
galaxy, the stellar population deduced is unusual but not extreme. 



6. Conclusions 

NGC 1222 is a remarkable galaxy. While there are uncertainties in the models we used 
for analyzing ionized gas emission, the results are clear and not model-dependent: NGC 
1222 has the highest excitation nebular spectrum and the highest deduced upper mass limit 
for the starburst stellar mass function of any galaxy yet studied that is not a metal-poor 
blue dwarf. NGC 1222 has about the same infrared diagnostic line ratios as the iconic 
metal-poor dwarf NGC 5253, but NGC 1222 has near-solar metallicity. It is unique among 
currently studied infrared galaxies. 

The infrared-line ratios of NGC 1222 would be extreme, for a galaxy of near-solar 
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metallicity, even without the extra factor added by our spatial modehng. With the 
models, the line ratios inferred are unique. That the inferred line ratios do resemble those 
observed in low-metal dwarfs suggests that perhaps the infrared emission is generated 
in a low-metallicity region. It would explain the results if, for example, NGC 1222 
is actually a lo w-metal content dwarf t hat h as merged with a system of higher metal 



content, and if 



Petrosian fc Burenkov 



(119931 ) measured their optical spectra in an area 



heavily contaminated with metals from the second system. The fact that there is a bright 
companion within 2 kpc of the central source of NGC 1222 indicates that the galaxy is 
currently experiencing a significant interaction. It would also would be consistent with 
another anomaly of NGC 1 222: its low molecular gas content (attempts to m easure the CO 



have only set upper limits: 



Chini. Krugel. fc Steppe 



19921: 



Elfhag et al. 



1996h . Starburst 



galaxies are usually rich in molecular gas unless they are metal-poor dwarfs. But NGC 1222 
cannot be a mis-classified dwarf galaxy, because it is too big and too bright. It has roughly 
half the infrared flux of NGC 5253, yet is 10 times further away. Its total luminosity puts it 
well outside the range of a dwarf galaxy. 

Higher spatial resolution measurements of the mid-infrared lines, specifically [Nell] 
and [SIV] which can be observed from the ground with sub-arcsecond resolution, would 
permit us to see directly the compact and extended components, and to determine how 
much extinction affects each. Higher resolution maps in the radio would find the spatial 
structure. The source cannot be a single cluster: at 1.3 cm its diameter is over 500 pc and 
its ionization is equivalent to 1.4 x 10^ 07 stars, twenty times that of the "supernebula" in 
NGC 5253. If the ionization is produced by stars with a Salpeter IMF and a lower mass 
cutoff of 1 M0 the total mass in stars is 2 — 3 x 10'' Mq. This exceeds by a large factor the 
most massive star clusters known. The high 1.3 cm flux raises the possibility that there are 
very compact, dense, rising-spectrum sources in NGC 1222. This would put it in the class of 
the most extreme and youngest starbursts, sources that are currently not well understood. 
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Why is this galaxy so different from most solar metallicity galaxies? The answer may 
be in the two companions which it is either absorbing or otherwise interacting with. A 
merger with one companion usually results in an intense burst of star formation. Perhaps 
in a double merger the process becomes especially extreme. 
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Fig. 1. — Radio maps of NGC 1222. top left: 6 cm image. Contour levels are 
0.18 mJy/beam. top right: 1.3 cm image. Contour levels are ±2"/2x 0.35 mJy/beam. The 
Spitzer IRS slit is superimposed, bottom: Contours of 1.3 cm emission on grey scale of 6-cm. 
Grey-scale is in mJy/beam, shown at top. Contour levels are ±2"/^x 0.15 mJy/beam. 
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Fig. 2.— The spectrum of NGC 1222 observed with IRS on SPITZER, reduced with SPICE. 
Order overlaps and noise spikes have been removed. 
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Fig. 3.— Line ratios calculated from STARBURST99+MAPPINGS simulations with ob- 
served line ratios. From top to bottom starting at the left of each figure are the derived 
ratios as a function of time for starbursts with upper mass limits of 100, 70, 50, 40, and 
30 Mq. Horizontal lines are the ratios observed in NGC 1222 for (top to bottom) Model 2, 
Model 3, Model 1. 




012345678 

Time after starburst (Myr) 

Fig. 4. — Line ratios found with CLOUDY, from Rigby & Rieke 2004. Horizontal lines are 
the observed NGC 1222 line ratios for Model 3 (top line) and Model 1. Model 2 (the most 
concentrated) would in all cases be off the top of the graph. 



-28- 



Table 1. Radio Observations 



Band 


Program 


Date 


[Beam Size 


r.m.s. 


Ghz 






// 


mjy 


C, 5 


AS286 


11/23/87 


1.75 X 1.2 


0.05 


K, 22 


AT309 


6/24/05 


0.96 X 0.53 


0.17 



Table 2. Infrared Lines in NGC 1222 



Observation 


[Nelll] 15.55//m 


[Nell] 12.81/im 


[SIV] 10.51//m 


[SIII] 18.71;im 




Jy 


Jy 


Jy 


Jy 


9071872-0006-5-E1394294 


1.7 


1.35 


0.32 


1.1 


9071872-0007-5-E1394295 


1.72 


1.6 


0.38 


1.37 



Table 3. Modeled Line Ratios: Compact Source 



Spatial Model [NeIII]/[NeII] [SIV] /[SIII] [SIV] /[Nell] 



Model 1 1.17 0.28 0.24 

Model 2 ^ 10.4 2.6 2.1 

Model 3 " 2.6 0.63 0.53 



^The galaxy is treated as one source 

''90% of [Nelll], [SIV] and 10% of [Nell] and [SIII] in compact 
source. 

^=90% of [Nelll], [SIV] and 40% of [Nell] and [SIII] in compact 
source. 
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Table 4. Line Ratios: Diffuse Emission 



Spatial Model [NeIII]/[NeII] [SIV]/[SIII] [SIV]/[NeII] 



Model 1 1.13 0.25 0.26 

Model 2 0.14 0.02 0.026 

Model 3 0.19 0.04 0.04 



